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We report first-principles calculations of the effects of quasiparticle self-energy and electron- hole 
interaction on the optical properties of single-walled BN nanotubes. Excitonic effects are shown 
to be even more important in BN nanotubes than in carbon nanotubes. Electron-hole interactions 
give rise to complexes of bright (and dark) excitons, which qualitatively alter the optical response. 
Excitons with binding energy larger than 2 eV are found in the (8,0) BN nanotubes. Moreover, 
unlike the carbon nanotubes, theory predicts that these exciton states are comprised of coherent 
supposition of transitions from several different subband pairs, giving rise to novel behaviors. 
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Boron nitride nanotubes (BNNTs) are isoelectronic 
to carbon nanotubes (CNTs); however, their electronic 
properties are quite different. Whereas carbon nan- 
otubes are metals or semiconductors with different size 
bandgaps depending on diameter and chirality Pj, BN 
nanotubes are wide gap insulators 0, 0. Although BN- 
NTs have been synthesized since 1995 only recently 
optical measurement on single-walled BNNTs has been 
performed Theoretical calculations 0, , as well as 
experiments 0, ^3 ? have shown that excitonic effects 
dramatically alter the behavior of the optical response of 
single- walled CNTs. For the BNNTs, these effects are ex- 
pected to be even more important due to the wide band 
gap nature of BNNTs. 

Experimentally it is found that BN nanotubes favor 
zigzag structure in current synthesis processes Thus, 
we focus our study on the zigzag tubes. Our calcula- 
tions on the (8,0) single- walled BNNT show that, indeed, 
many-electron effects lead to the formation of strongly 
bound excitons of multi-band character with extraordi- 
narily large binding energies, which dramatically change 
its optical absorption spectrum. 

To compute the optical response, we use the method of 
Rohlfing and Louie jl^ in which electron-hole excitations 
and optical spectra are calculated from first principles in 
three steps. First, we treat the electronic ground state 
with ab initio pseudopotential density-functional theory 
(DFT) pj. Second, we obtain the quasiparticle energies 
^nk within the G W approximation for the electron self- 
energy S ^] by solving the Dyson equation: 

+ Vion + VHartree + 5](£^nk) 
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Finally, we calculate the coupled electron-hole excita- 
tion energies and optical spectrum by solving the Bethe- 
Salpeter (BS) equation of the two-particle Green's func- 
tion lillii: 

k'^''c' 

where A^^^ is the exciton amplitude, K^^ is the electron- 
hole interaction kernel, and |ck) and \vk.) are the quasi- 
electron and quasihole states, respectively. 



The DFT eigenvalues and wave functions were ob- 
tained within the local density approximation (LDA) 
using a plane-wave basis 16] with an energy cutoff of 100 
Ry. ab initio Troullier-Martins pseudopotentials [l^ in 
the Kleinman-Bylander form jJJ] were used (with cutoff 
Tc for boron and nitrogen of 0.79 A and 0.63 A, respec- 
tively). For convergent results to better than 0.05 eV, 
up to 32 k points in the one-dimensional Brillouin zone 
were used for the GW calculations and for solving the 
BS equation. All calculations were carried out in a su- 
percell geometry with a wall-to-wall intertube separation 
of 9.5 A to mimic isolated tubes, together with a trun- 
cated Coulomb interaction in the radial direction in order 
to eliminate unphysical interactions between periodic im- 
ages on the different tubes. As shown in Ref. it is 
important to truncate the Coulomb interaction because 
these unphysical interactions would increase the effective 
screening in the system and hence reduce both the self- 
energy correction and the exciton binding energy. Be- 
cause of depolarization effects in nanotubes [2^, strong 
optical response is only observed for light polarized along 
the tube axis {z). We consider here only this polariza- 
tion. 

Figure ^a) shows the quasiparticle energy corrections 
to the LDA energy eigenvalues. First, we note that these 
corrections are quite large, in comparison to those for 
bulk hexagonal BN (h-BN) and SWCNTs. The quasi- 
particle corrections open the LDA gap of bulk h-BN by 
~1.58 eV near zone center or the F-point [i^l, while the 
gap opening in the (8,0) SWBNNT near the F-point is 
~3.25 eV. This is a consequence of enhanced Coulomb 
interaction in reduced dimension 7]. Also, due to its 
larger gap, which weakens screening, the quasiparticle 
corrections to the gap in the (8,0) SWBNNT are larger 
than those for a similar SWCNT (which are ~ 1.15 eV 
near the F-point 0). Second, the quasiparticle correc- 
tions have a complex band- and energy-dependence, so 
for accurate results they cannot be obtained by a simple 
scissor shift operation. The corrections depend on the 
character of the wavefunction. For example, states of the 
fourth lowest conduction band in the LDA bandstructure 
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FIG. 1: Difference between the CVF quasiparticle energy and 
the LDA Kohn-Sham eigenvalue plotted as a function of the 
energy of the states (a) and quasiparticle band-structure (b) 
for the (8,0) SWBNNT. Empty circles in (a) and the dashed 
line in (b) show the nearly- free-electron tubule states. 
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FIG. 2: (color online). Absorption spectra of the (8,0) 
SWBNNTs. The imaginary part of the polarizability per tube 
0^2 (cj) is given in unit of nm^. (See text.) The spectra are 
broadened with a Gaussian of 0.0125 eV. 



are localized inside the tube, and are nearly free-electron- 
like states. These tubule states form a separate branch 
in the quasiparticle correction diagram with significantly 
smaller corrections. Figure^b) depicts the quasiparticle 
bandstructure of the (8,0) SWBNNT. The arrows indi- 
cate the optically allowed interband transitions between 
four pairs of bands which give rise to the lowest-energy 
peak structures in the non-interacting optical spectrum 
in Fig. 13 (The fourth lowest energy conduction band in 
the LDA bandstructure becomes the lowest energy con- 
duction band in the quasiparticle bandstructure.) 

FigureI21depicts the optical absorption spectrum calcu- 
lated with and without electron-hole interaction effects. 
The plotted quantity is the imaginary part of the calcu- 
lated dielectric susceptibility, X = (^ ~ l)/47r, multiplied 
by the cross-sectional area of the supercell perpendicu- 
lar to the tube axis. This quantity a, as defined above, 
gives the polarizability per single tube in units of nm^; 
so the susceptibility of an experimental sample contain- 
ing a density of n infinitely long tubes per unit area may 
be obtained as x = ^cif- The absorption profile changes 
dramatically when the electron-hole interaction is taken 
into account. We use the label /, and // to denote 
the excitons that give rise to the absorption peaks in 
the figure. Subscripts 1, 2, 3, and 4 refer to the ground, 
first-excited, second-excited, and third-excited states of a 
particular bright exciton series, respectively. We observe 
a strong absorption peak at 5.72 eV, which corresponds 
to a bound exciton (/i) with a binding energy of 2.3 eV. 
The area under this peak, which may be used in compar- 
ing with experiment, is 0.87 nm^eV. Excitons Ii and /( 
are different states, made up of transitions from the same 
set of four pairs of valence and conduction subbands of 
the (8,0) BNNT, all of which have similar quasiparticle 
transition energies from 8.1 eV to 8.3 eV (See arrows 
in Fig. ^b)). These transitions are coupled strongly to 
each other by the electron-hole interaction to form the 



lowest optically active states (the singly-degenerate /i 
and doubly-degenerate /{). This behavior is very differ- 
ent from the (8,0) SWCNT in which the exciton states 
are composed mainly of transitions between a single pair 
of quasiparticle bands. 

The mixing of transitions of different subbands alters 
the electron-hole wavefunction - localizing further the 
electron amplitude with respect to the hole position in 
real space and making it deviate from a 1-dimensional be- 
havior with spatial variations in directions perpendicular 
to the tube axis. Figure Ol^a) shows the isosurface plots 
of the electron distribution |<l>(re, r/j,)p with the hole po- 
sition Yh fixed (the black star in the figure) for the first 
bound exciton (/i). Figure Ofb) quantifies the electron- 
hole correlation for this state by plotting \^\^ along the 
tube axis after integrating out the electron coordinates 
in the perpendicular plane (the hole position is set at 
zero). The position of the peaks in FiglH^b) corresponds 
to the position of plane of boron atoms, i.e., the pho- 
toexcited electron is localized on the boron atoms near 
the hole. Thus, as expected, the photo-excitation pro- 
cess corresponds to a transfer of electron from nitrogen 
atoms to nearby boron atoms; but the resulting electron 
and hole amplitudes are strongly correlated with an ex- 
tent of only a few inter-atomic distances. Figure[2fc) is a 
cross-sectional plot, showing the excited electron proba- 
bility distribution in a plane that is perpendicular to the 
tube axis and contains the hole as well as other nitrogen 
atoms. 

As a comparison to carbon nanotubes. Fig. intd)-(f) 
shows similar quantities as in Fig. inta)-(c) but for the 
first bright bound exciton in the (8,0) SWCNT 0. In 
the figure, the hole is fixed slightly above a carbon atom. 
The exciton in the (8,0) SWBNNT is significantly more 
tightly bound than that in the (8,0) SWCNT and cannot 
really be viewed as a 1-D object. The root-mean-square 
size of the exciton along the tube axis is 3.67 A for the 
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FIG. 3: (color online), (a)-(c): Wavefunction of the low- 
est energy bright exciton of the (8,0) SWBNNT. (a) Isosur- 
face plot of electron probability distribution |$(re,r^)|^ with 
the hole fixed at the position indicated by black star, (b) 
|$(re,r^)|^ averaged over tube cross section. Hole position 
is set at zero, (c) |$(re,rh,)|^ evaluated on a cross-sectional 
plane of the tube, (d)-(f): Wavefunction of the lowest energy 
bright exciton of the (8,0) SWCNT. Plotted quantities are 
similar to those in (a)-(c). 



(8,0) SWBNNT and 8.59 A for the (8,0) SWCNT, and 
their binding energies are 2.3 eV and 1.0 eV, respectively. 
This difference in behavior is due to the wide bandgap 
and weaker screening in SWBNNT. Also, we note that 
while the binding energ y of the excitons in the bulk h- 
BN is only 0.7 eV [23,131, the binding energy in the (8,0) 
SWBNNT is more than three times larger. 

Figures 4(a) and 4(b) show similar quantities as in 
Fig. ISl^b) for the excitons /{ and /2- For exciton /(, 
the electron is less tightly bound to the hole than in 
exciton Ii. The state /2, which is an excited state of 
exciton /i, is also more diffuse than Ii and the electron 
amplitude is not at a maximum near the hole which is 
the case for /i(Fig. Oljb)). We also note that, for the 
(8,0) SWBNNT, there are numerous dark excitons dis- 
tributed rather uniformly in energy below and among the 
bright excitons shown in Fig. [21 The energy of the lowest 
doubly-degenerate bound dark exciton is at 4.63 eV. This 
dark exciton is made up of transitions from the highest 
valence band to the lowest conduction band (the NFE 
tubule state) in the quasiparticle bandstructure, and has 
a binding energy of 1.94 eV with respect to these inter- 
band transition energies. 

The various lowest-energy exciton states (for both 
bright and dark excitons) derived from the various dif- 
ferent sets of interband transitions, on the average, have 
a large binding energy of about 1.9 eV. However, the 
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FIG. 4: (color online). Wavefunctions of excitons of the 
(8,0) SWBNNT. Plotted quantities are similar to those in 
Fig. 3(b). 



binding energy of the first bright exciton is 2.3 eV. We 
ascribe this extra binding energy of about 0.4 eV to the 
fact that four different sets of interband transitions are 
strongly coupled in forming the first bright exciton Ii. 
This strong coupling mixes states from the different tran- 
sitions, splits the excitation energy levels, and increases 
the binding energy of the final lowest-energy exciton. 

Lauret et. al. 5] recently measured the optical proper- 
ties of samples containing SWBNNTs and observed three 
absorption peaks at 4.45, 5.5, and 6.15 eV, respectively. 
These authors assigned the 4.45 eV peak to originate 
from a bound exciton state. In our calculation, the first 
strong absorption peak is at 5.72 eV which is about 1.3 
eV blue shifted from the observed peak. A direct compar- 
ison between theory and experiment however should not 
be made. The average diameter of the tubes measured 
in the experiment is 1.4 nm while the diameter of the 
relaxed (8,0) SWBNNT in our work is 0.65 nm. Also, 
in the experiment, there would be tubes with different 
chiralities. 

It is unlikely that the large discrepancy of 1.3 eV is due 
to differences in diameter or chirality. At the LDA level, 
which neglects both self-energy corrections and excitonic 
effects, the position of the first broad absorption peak 
for SWBNNTs with diameter of 0.48 nm for a (6,0) tube 
and 2.13 nm for a (27,0) tube is at 4.8 eV and 5.6 eV, 
respectively |2^. The difference between them is only 
0.8 eV. Furthermore, assuming same chirality, larger di- 
ameter tubes tend to have their first absorption peak at 
higher energy. For these reasons, it is difficult to explain 
the 1.3 eV discrepancy. The calculated first peak po- 
sition for the (8,0) tube is however rather close to the 
observed 5.5 eV peak. Also, the observed difference be- 
tween the second and the third absorption peak position 
in the experiment is 0.65 eV, very close to the difference 
between the first and the second absorption peaks in our 
calculation which is 0.62 eV. But the difference between 
the first and the second peak position in the experiment 
is 1.05 eV. Another difference is that the theory is for 
an isolated tube, while experimentally the tubes are sur- 
rounded by a dielectric medium which can modify the 
excitation energies. We thus suspect that the observed 
second peak at 5.5 eV may be due to an exciton, corre- 
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spending in nature to our first absorption peak at 5.72 
eV. However, an additional complication is the existence 
of the many dark excitons, which could be made active 
by external perturbations. The exact interpretation for 
the experimental observations can therefore only be made 
after further experimental and theoretical studies. 

For SWCNTs, the effect of the surrounding dielectric 
medium on the optical spectrum is expected to be small 
(even though it can be important for the exciton bind- 
ing energy) due to an almost cancellation between the 
quasiparticle self-energy correction and the binding en- 
ergy of excitons 7]. For SWBNNTs, screening by ex- 
ternal medium would be more important because intrin- 
sic screening is much weaker. In particular, in isolated 
SWBNNTs, we find that the energy difference between 
the quasiparticle self-energy correction and the binding 
energy of the exciton is large (about 0.9 eV). In the pres- 
ence of a dielectric medium, we expect this energy differ- 
ence to decrease, which would result in a decrease in the 
excitation energies. 

Among previous theoretical works on the optical prop- 
erties of BNNTs, Guo and Lin carried out calcula- 
tions at the LDA-RPA level without considering many- 
electron effects. Their optical absorption spectra are 
qualitatively very different from the present final results. 
From their results for the (6,0) and (9,0) tubes (showing 
broad absorption peaks at 4.8 eV and 5.0 eV, respec- 
tively), we can deduce a peak position for the (8,0) tube 
to be near 4.9 eV, as we find in our LDA-RPA level calcu- 
lation. The first peak position in Fig. [21 with electron- hole 
interaction included is blue shifted by about 0.9 eV from 
that of the LDA-RPA calculation. Ng and Zhang on 
the other hand, used a time-dependent, localized-density- 
matrix approach based on a semiempirical Hamiltonian 
to compute the optical properties of SWBNNTs. In their 
work, the absorption spectrum of (8,0) BNNT has a very 
broad first peak at an even higher energy than ours, at 
near 6.2 eV. 

In summary, we have done calculation on the zigzag 
(8,0) SWBNNT to study the effects of many-electron in- 
teractions on its optical response. The Gl^ corrections 
to the quasiparticle excitation energies of the SWBN- 
NTs are significantly larger than those for SWCNTs or 
bulk h-BN. Also, the quasiparticle energy corrections are 
found to be complicated so that interpolation by a sim- 
ple scissor shift operation is not a good scheme for accu- 
rate calculation. Theory predicts that, unlike the non- 
interacting case, the absorption spectrum of the (8,0) 
SWBNNT is dominated by a huge peak at 5.72 eV, due 
to an exciton with a large binding energy of 2.3 eV. This 
exciton state is made up of optically-allowed transitions 
between four different pairs of subbands. Moreover, an 
intricate set of dark excitons is found to exist. Self-energy 
and electron-hole interaction effects therefore are even 



more important in the optical response of the SWBN- 
NTs than in the SWCNTs. 
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